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Environments of z > 5 quasars: searching for protoclusters at 
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ABSTRACT 

The most massive haloes at high redshift are expected, according to hierarchical cosmolo- 
gies, to reside in the most biased density fields. If powerful active galactic nuclei (AGN) are 
expected to exist anywhere in the early Universe (z > 5), it is within these massive haloes. 
The most luminous of these AGN, powered by supermassive black holes (SMBHs) ~10 9 M Q , 
thereby present an opportunity to test models of galaxy formation. Here, we present submil- 
limetre (submm) continuum images of the fields of three luminous quasars at z > 5, obtained 
at 850 and 450 /im using the Submm Common-User Bolometer Array (SCUBA) on the James 
Clerk Maxwell Telescope (JCMT). N-body simulations predict that such quasars evolve to be- 
come the central dominant galaxies of massive clusters at z = 0, but at z = 5 — 6 they are 
actively forming stars and surrounded by a rich proto-filamentary structure of young galaxies. 
Our purpose in taking these images was to search for other luminous, star-forming galaxies 
in the vicinity of the signpost AGN and thus associated with such a protocluster. Two of the 
quasar host galaxies are luminous SMGs in their own right, implying star-formation rates 
(SFRs) ^10 3 M yr _1 . Despite the coarse 850-/im beam of the JCMT, our images show evi- 
dence of extended emission on a scale of ^100-kpc from at least one quasar - indicative of a 
partially resolved merger or a colossal host galaxy. In addition, at >3a significance we detect 
12 (5) submm galaxies (SMGs) at 850 /im (450 /im) in the surrounding fields. Number counts 
of these SMGs are comparable with those detected in the fields of z ~ 4 radio galaxies, and 
both samples are systematically overabundant relative to blank-field submm surveys. Whilst 
the redshift-sensitive 850 /im/450 /im and 850/im/1.4GHz flux density ratios indicate that 
some of these SMGs are likely foreground objects, the counts suggest that many probably lie 
in the same large-scale structures as the quasars. 
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1 INTRODUCTION 

The search for z > 5 galaxies ha s taken on a new urgency in recent 
years. The epochal discovery of iGunn & Petersori 1 19651) absorp- 
tion in the spectra of z > 6 quasars teecker et all200lh suggests 
that reionisation was complete by z — 6. On the other hand, po- 
larisation results from the Wilkinson Microwave Ani sotropy Probe 
( WMAP) point to the onset to reionisation at z ~ 12 ( Sperg eTet al.l 
200-J. A primary goal of observational cosmolog y is, conse- 
quently, to locate and to cha racterise the galaxies (e.g. Bunker et al.l 
l2004lBouwens et al.120061) . quasars (e.g. lFan et alj2006l) and - via 
gamma-ray bursts - even individual stars ( Jako bsson et al. 2006 ; 
lHaislip et al]|2006t) that existed during this period of cosmic phase 
transition, ~0.4- 1 Gyr after the Big Bang. 

Although it is thought that the ultraviolet (UV) light from 
young stars, rather than AGN, provided the bulk of the ionizing 



radiation dFan et al 1 l2006l : ISrbinovskv & Wvi1faell2007t) . very high- 
redshift quasars are nevertheless important objects in their own 
right. SMBHs are found ubiquitously in the cores of nearby galax- 
ies; studying their accretion-powered growth phase, during which 
they manifest themselves as high-redshift AGN, provides a win- 
dow on the formation and evolution of the ancestors of present- 
day massive galaxies. The very existence of luminous AGN at 
redshifts above five imposes important constraints on processes 
governing the formation of massive galaxies, SMBHs and proto- 
clusters in the early Universe. Their bolometric powers imply ac- 
cretion onto SMBHs at a prodigious rate (~lMoyr _1 , assum- 
ing Eddington-limited accretion). Measurements base d on the line 
width of Mg II indicate SMBH masses of ~10 9 5 M jWillott et all 
120031). implying a host dark matter halo of mass 
i Croometal]|2005h . Submm photometry of the host galaxies of 
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Table 1. Properties of target z > 5 quasars and their host galaxies. 



Source Name 


2 


M B 


R.A. Dec. 
(J2000) 


ophot 
'-'l250Mm 

(mJy) 


crphot 

(mJy) 


Qphot 
' 5 450nm 

(mJy) 


•51.4GHz 

d«Jy) 


SDSS J0756+4104 
SDSS J1030+0524 
SDSS J1044-0125 


5.09 
6.28 
5.73 


-26.9 
-27.7 
-28.0 


07 56 18.14 +41 04 08.6 
10 30 27.10 +05 24 55.0 
10 44 33.04 -01 25 02.1 


5.5 ±0.5 

<3.4 
2.5 ±0.55 


13.3 ± 2.1 
1.3 ±1.0 
6.1 ±1.2 


14±19 
-21 ±10 


65 ±15 
<61 
<79 



z > 5 quasars reveal copious quantities of dust (~ 10 s 9 Mq - 
IPriddev et ai]|2003bl ; iRobson et aill2004h and CO emission line s 
imply ~ 10 10 " 11 Mq of molecular gas dWalter et alfeOOl 120041) . 
Both are indicative of strong star formation; indeed, a colossal 
starburst is required if the dust has been synthesised on the short 
timescale available (<1 Gyr since the Big Bang; <600Myr since 
the onset of reionization). 

The existence of massive, collapsed haloes at high redshift 
might appear to conflict with the prevailing hierarchical cold dark 
matter (CDM) model of structure formation, in which the initial 
density fluctuations have higher power on smalle r scales. Using the 
Press- Schechter approximation, for example, lEfstathiou & Reesl 
(1988) predicted a marked decline in the number density of 
quasars at z > 5. To accommodate the growing sample of lu- 
minous AGN spectroscopically confirmed to lie at z > 5 within 
the CDM scheme, they are identified as corresponding to rare, 
high-cr peaks in the initial overdensity distribution. As such, they 
constitute excellent laboratories in which to test the predictions 
of structure-formation models operating in extremis. Specifically, 
high-a haloes should be biased - more spatially correlated than 
the underlying mass - since neighbouring fluctuations sitting atop 
a low-amplitude, large-scale mode of the ov erdensity spe ctrum 
are boosted over the collapse threshold sooner dKaiserlll984h . The 
vicinity of a high-redshift, luminous AGN would thus be expected 
to contain an overabundance of massive, actively star-forming 
galaxies, the progenitor of a rich cluster at z = 0. More recently, 
the phenomenon of "downsizing" provides a framework in which 
galaxy formation at high redshift is completed earlier in more 
massive objects. Specificall y, the semi-analytic "Anti-hierarchical 
Baryon Collapse" model of [Granato et al.l f2004) suggests an inti- 
mate evolutionary link between quasars and SMGs, feedback from 
the AGN playing an important role in regulating star formation in 
the host. 

The submm waveband is an eminently appropriate place in 
which to search for structure, as traced by dusty, star-forming 
galaxies surrounding high-redshift AGN. This conjecture is sup- 
ported by submm/mm imaging campaigns which have identified 
an excess of SMGs in fields centred on z ~ 4 radio galax- 
ies dlvison et al.ll2000al ; Istevens et ai1l2003l ; be Breuck et al J2004l; 
Greve et alj 120071) and a z ~ 2 absorbed QSO dStevens et al.l 
2004h . n this paper, we extend this work to higher reds hift, and ex- 
pand upon our earlier photo metric submm programmes dlsaak et all 
l2002l ; |Priddev" et al. 2003a), obtaining very deep submm imaging of 
the fields of three z > 5 radio-quiet quasars. 

Throughout, we assume cosmological param eters Qm = 
0.27, Q.a = 0.73 and H = 71kms" 1 Mpc" 1 dSpergel et all 
120071) . 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 The sample 

Three fields centred on quasars at z > 5 were selected for obser- 
vation. The quasars we re all discovered as par t of t he Sloan Digital 
Sky S urvey (SDSS) bv lFan et al.l d2000l.l200lh and lAnderson et al.l 
Their properties are listed in TableQ] all three are luminous 
at optical wavelengths (absolute B magnitudes in Table Q] are ex- 
trapolated from 1450A continuum flux assuming a spectral index, 
—0.5), implying black hole masses of ~1Q 9 3 ~ 9 7 Mq, assuming a 
bolometric correction from the B band of 12 dElvis et "all ll994l) 
and Eddington-limited accretion. The sample was also designed 
to provide contrast in f ar-infrared luminosity, containing both the 
brightest member of the lPriddev et al.l d2003bh submm photometry 
sample (SDSS J0756+4104, Sgso^m = 13.4 mJy) and the deepest 
upper limit (SDSS J1030+0524, Ssso^m < 2mJy). 

2.2 Observations 

Jiggle-m ap observations wer e made using the SCUBA bolome- 
ter array dHolland et alJI 1999h on the 15-m JCMT over the period 
between 2001 Feb 28 and 2002 Mar 23. In each case the array 
was centred on the optical position of the quasar. Every 16 s the 
telescope was nodded in either right ascension or declination by 
30 arcsec in an ON-OFF-OFF-ON sequence, with the secondary 
mirror chopping continually at ~7 Hz between the ON and OFF po- 
sitions. For SDSS J1044— 0125, we nodded and chopped east-west 
(E-W) throughout. For SDSS J0756+4104 and SDSS J1030+0524 
we nodded and chopped north-south (N-S) as the target rose and 
set, and E-W when the target was near transit. By adopting this lat- 
ter strategy we ensured the chop direction was as close to azimuthal 
as possible and that the region within 30 arcsec of the outer bolome- 
ters was chopped onto the array; also, using several position angles 
reduced the chance of chopping one submm source systematically 
onto another. 

Dynamic scheduling at JCMT ensured weather conditions 
were exceptionally good for our observations. The sky was stable 
and very transparent, with a median t^sghz °f around 0.043 (often 
<0.04, never worse than 0.09). Flux calibration was determined by 
imaging CRL 618, Mars and Uranus and is accurate to 10 percent at 
850 /im and 25 percent at 450 /jm. A long imaging scan of 3C 345 
was obtained, using the same observing strategy, to establish an 
accurate point spread function (PSF). 

The total number of integrations/time expended on each 
field were: SDSS J0756+4104 (304 integrations; 38.9 ks); SDSS 
J1030+0524 (280 integrations; 35.8 ks); SDSS J1044-0125 (312 
integrations; 39.9 ks). 

In addition, we publish here for the first time a 250GHz on- 
off photometric observation of SDSS J1044— 0125 obtained using 
the Max Planck Institut fur Radioastronomie Millimetre Bolome- 
ter (MAMBO) array on the Institute de Radioastronomie Mil- 
limetrique (IRAM) 30-m telescope on Pico Veleta, Spain. The sky 
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opacity during the observations was in the range 0. 164—0.260. Data 
were reduced using the GILDAS software package, yielding a 1.2- 
mm flux density 2.5 ±0.6 mjy (4.5a). 



2.3 Data reduction 

Data from SCU BA were reduced in a sta ndard manner with the 
SURF package dJenness & Lightfoot|[l998l) . using modifications to 
the SETBOLWT and REBIN tasks developed for the s urvey of high- 
redshi ft radio galaxies by_ Ivison et all d2000al) and IStevens et al.l 
(2003), as described bv llvisonl d2006h . Briefly, we created an astro- 
metric grid of 1-arcsec 2 square pixels then determined accurately 
weighted and calibrated signal and noise measurements for each 
pixel, each measurement being wholly independent of values for its 
neighbouring pixels, with 13.4-arcsec FWHM resolution at 850 /im. 
These signal and noise values reflect the stream of data collected 
when bolometers are centred within the region of sky correspond- 
ing to a particular pixel. These unsmoothed images were used for 
source detection (see ij2.4b . 

The submm images of each of the three quasar fields - 850-/im 
contours superposed on greyscale representations of the 450-/im 
data - are shown in Fig.Q] where we have smoothed with 6- and 3- 
arcsec FWHM Gaussians and the resulting 850- and 450-/xm maps 
have typical noise levels of 1.5 and 7 mjy beam -1 . For comparison, 
the 5a 850-/im confusion limit for the JCMT is m 1-1.5 mjy, so we 
have reached close to the depth at which confusion is likely to begin 
posing problems. 



2.4 Source detection 

Source detection at 85 /tm was accom plished using the algorithm 
described in detail by IScott et ail d2002h . utilising the signal and 
noise maps and a PSF (measured for the blazar, 3C 345) to perform 
a simultaneous maximum-likelihood fit to all the potentially signif- 
icant peaks in each map. Sources down to a significance level of 
3 a are listed in Table [2] 

Adopting a 3a detection threshold, at 850 pm (450 /im) a total 
of 14 (7) sources are detected: 2 (1) QSOs and 12 (6) companions in 
the three fields, with flux densities between 4.4 (16) and 13.4 mjy 
(56mJy). Five sources, including one QSO, are detected at both 
wavelengths (Table |2j- 

Discussing each of the fields in turn: 

SDSS J0756+4104 (z = 5.09) was detected via pho tometry-mode 
(on-of f) observations at 850 pirn: 13. 4 ±2.1 m jy ( Pridde v et al.l 
l2003bh and at 1200 /an: 5.5 ± 0.5 mjy dPetric et al.ll2003h . Submm 
emission is clearly visible in the maps (Fig. QJ at the opti- 
cal coordinates of the quasar (the map centre), at both 850 
(11.2 ± 1.0 mjy) and 450 pm (16±5mJy). To the eye, the source 
appears elongated; the best-fit 2-D Gaussian corresponds to 
(21.4 ±9.0) arcsec x (13.6 ±3.8) arcsec at a position angle (PA) 
of 69 ±37°. 

Three further sources were detected (>3a) in the SDSS 
J0756+4104 field at 850 /im, each close the edge of the map (one 
beyond the areal coverage at 450 /im). In addition, two sources 
were detected at 450 /im, though neither has a robust 850-/im coun- 
terpart (3a < 4.5 mjy). This casts some doubt on their reality; cer- 
tainly, if they are real emitters, they are likely to be at a much lower 
redshift than the central quasar. 

SDSS J1030±0524 (z = 6.28) has been the object of consider- 
able observational study. IStiavelli et alj d2005l) . for example, de- 
tected two candidate z ~ 6 sources in Hubble Space Telecope 
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Figure 1. submm images of fields centred on three optically selected, radio- 
quiet quasars at z > 5. The 850-/im maps are plotted as contours, at levels 
-3, +3, 4, 5, 6... Xcr. The 450-/tm images are displayed as greyscales. 850- 
fim sources (at ^3<r significance) are identified with open crosses; their 
labels correspond to those in column 2 of Tablefl] Sources detected only at 
450 fim are circled. 
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Table 2. Properties of SMGs detected in quasar fields. 



IAU Source Name 


Id in 




R. 


A. 


Dec. 




S/N 


5450/jm 




Fig.E 






(J2000) 






±a (mjy)t 




±a (mjy)t 


SDSS J0756+4104 


#1 


07 


56 


18.15 


+41 


04 


07.5 


13.4 ± 1.0 


13.9 


16±5 


olvllVl J U 1 JO 1 J .Jjt4 1 UjUo. / 


#2 


07 


JO 


15.55 


+41 


03 


08.7 


A A 4- 1 A 


^ [ 


3<r < 24 


kJlVllVl J U / JU 1 J .OUTt 1 U JUU.i 


#3 


07 


56 


15.86 


+41 


03 


00.2 


4 6 ± 1 4 




3cr < 24 


kjlVllvl JU / JO 1 7.U7Tt 1 UJJH-.4 


#4 


07 


56 


19.09 


+41 


05 


34.4 


13 1 ± 3 4 


T Q 


t 
+ 


ijlVllvl JU / JOZ 1 .oj+t 1 Lf+J? .y 


#5 


07 


56 


21.83 


+41 


04 


59.9 






^f. -4- 1 ^ 

JO ZL 1 J 


SMM J075623.58+410424.6 


#6 


07 


56 


23.58 


+41 


04 


24.6 


3cr < 5.1 


- 


25 ±8 


SMM J 103029.03+052502.9 


#1 


10 


30 


29.03 


+05 


25 


02.9 


7.5 ± 1.7 


4.3 


42 ±12 


SMM J103028.21+052452.3 


#2 


10 


30 


28.21 


+05 


24 


52.3 


5.3 ± 1.7 


3.3 


3cr < 33 


SMM J103025.41+052431.7 


#3 


10 


30 


25.41 


+05 


24 


31.7 


6.8 ± 1.7 


4.3 


49±11 


SMM J103022.04+052435.0 


#4 


10 


30 


22.04 


+05 


24 


35.0 


7.0 ±2.2 


3.3 


3cr < 48 


SMM J103026.29+052524.9 


#5 


10 


30 


26.29 


+05 


25 


24.9 


8.5 ± 1.7 


5.0 


45 ± 10 


SDSS J1044-0125 


#1 


10 


44 


33.15 


-01 


25 


03.6 


5.6± 1.0 


6.2 


3cr < 24 


SMMJ104433.31-012527.7 


#2 


10 


44 


33.31 


-01 


25 


27.7 


7.5 ± 1.5 


5.5 


24 ±7 


SMM J104429.22-012529.0 


#3 


10 


44 


29.22 


-01 


25 


29.0 


6.1 ±1.7 


3.5 


3cr < 24 


SMM J104430.74-012528.2 


#4 


10 


44 


30.74 


-01 


25 


28.2 


7.0± 1.7 


4.5 


3cr < 24 


SMMJ104430.71-012357.5 


#5 


10 


44 


30.71 


-01 


23 


57.5 


6.8±1.9 


3.0 


t 



f Errors exclude the uncertainty in absolute flux calibration: 10 (25) per cent at 850 fim (450 ^m). 
| Outside region covered by 450-/xm data. 



(HST)-ACS images of the field, which they speculate may belong 
to the same dark matter halo. The quasar itself appears similar in 
all respects to low-redshift AGN: submm photometry failed to de- 
tect the quasar to a sensitive limit of 3<j < 3.0 mjy at 850 /^m 
dPriddev et alJ2003bl ) and a deep X-ray spectrum taken with XMM- 
Newton dFarrah et al]|2004) finds a photon index V ~ 2.1 with no 



evidence for intrinsic absorption (Nu < 8 x 10 cm ). The lack 
of gas and dust, plus the near-canonical X-ray spectral slope, sug- 
gest that this SMBH and its host have already passed through their 
major formation phase, only ~900 Myr after the Big Bang. 
Our submm images (Fig.QJ confirm the submm quiescence of the 
quasar/host galaxy, with 3-cr limits of 5.1 and 30mJy at 850 and 
450 /im, respectively. Indeed, there is an interesting negative fea- 
ture just south-west of the quasar corresponding to the western 
"off" position of the weakest SMG in this field, probably the pri- 
mary reason SMM J103028. 21+052452. 3 achieves an overall sig- 
nificance Of ^3(7. 

The surrounding field is crowded: five sources at 850 ^im and three 
at 450 /im. Each of the 450-^im sources lies within a fraction of a 
beamwidth of significant 850-/im emission, confirming the reality 
of both. 

Is it surprising that the least active quasar host should have the 
most active environment? It is possible that, by targeting luminous, 
opfica/Zj-selected quasars, we are selecting very massive objects 
that are already close to the end of their star-forming phase. By the 
rules of "downsizing", the evolution of the most massive objects is 
completed earliest, so it should not be surprising to see compan- 
ions to fully-formed, high-redshift quasars in an earlier evolution- 
ary stage (see also Section|33]l. 

A 1.4-GHz V ery Large Array im age (a ~ 20 /iJy) of the field was 
obtained by dPetric et al. I l2003h . comparison with which enables 
us to pinpoint counterparts at ot her wavelengths with sub-arcsec 
precision (e.g. Ilvison et"al] fl998). In addition, constraints can be 
placed on the redshifts of those sources with radio counterparts, 
or meaningful limits. SMM J103029.03+0452502.9 lies 6.8 arcsec 
from a radio source at a = 10 h 30 m 26. s 0, 5 = ±05°25'19."69 
with an integrated 1.4-GHz flux density of 427 ± 52 /iJy (Andreea 



Petric, private communication), giving a spectral index, «i 5 4ghz 
= 0.52 ±0.04. If the submm and radio sources are th e same (and 
this i s statistically very probable - P < 0.05 - see Ilvison et al.l 
2002), the FIR-radio correlation for star -forming galaxies pre- 
dicts z = l.Qtol dCarilli & Yunl l200d : see Figure SMM 
J103025. 41+052431. 7 is near-coincident (separation, 4.4 arcsec) 
with the core of a limb-brigh tened Fanaroff-Riley II (FRll) radio 
galaxy (Fig. 2 of IPetric et al.l) - another statistically significant as- 
sociation, i.e. P < 0.05. This raises the strong suspicion, allied 
with their relatively high 450-/im flux densities, that some or all of 
the other nearby SMGs might lie in the foreground, associated with 
this radio galaxy, rather than at the redshift of the quasa r. Such lu 
minous rad io galaxies are commonly found in clusters (Best et al 
2003, 2007) and s ome clusters are known to host SMGs ( Best 200: 
Webb et alj|2005l) . 

SDSS J1044 0125 (z = 5.73) has an optical spectrum with 
a Civ absorption trough sugg esting that this i s a broad ab- 
sorpti on line (BAL) quasar jMaiolino et al.ll200ll ; iGoodrich et al 



LI- 

al. 

2; 



200 1 ) . Correspondingly, the quasar is weak in X-rays l Brandt et al 



ys 
G 



2001), consistent with the BAL interpretation (e.g. Gallag her et al 
2006). The quasar/host was detected at 850 and 1200 /im using 
on-off measurements with SCUBA and MAMBO respectively: 
6. 1 ± 1.2 mjy and 2.5 ± 0.6 mjy (see Section l2~2l . 
In our SCUBA image, the quasar is detected robustly at 850 
with a flux density of 5.7 ± 1.5 mJy, with weak evidence of ex- 
tended emission (Section l3.4t . Four further sources are detected 
at 850 £tm; one at 450/xm, corresponding to the brightest 850-/im 
emitter. 



3 DISCUSSION AND ANALYSIS 
3.1 Photometric redshifts 

We can place crude constraints on the redshift of any 850-/im emit- 
ters in our quasar fields using detections or limits at 1.4 GHz and/or 
450 /j,m. Much of the uncertainty in this proc edure rests in the in - 
trinsic diversity of SEDs exhibited by SMGs ( Ilvison etal . 2000b). 
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2 4 6 8 

Redshift 

Figure 2. Top: submm (850 (im) to radio (1.4 GHz) spectral index as 
a redshift indicator. The curves show em pirical cal i bration s from (thick 
solid, with zfclcr sho wn as thin solid) Carilli & Yun 12000) and (dashed) 
Also for comparison are plotted (sma ll points ) SMG s 
with spectroscopic redshifts and radio fluxes from Chapma n et alj j2005h . 
Sources in the SDSS J0756+4104 field are shown as triangles; in the SDSS 
J1030+0524 field, squares; in the SDSS J1044-0125 field, circles. Radio 
detections have filled symbols; upper limits are unfilled. Two of the quasars 
(SDSS J0756+4104 and SDSS J1044-0125) are plotted at their spectro- 
scopic redshifts. A tentative radio counterpart to an SMG in the SDSS 
J 1030+0524 field (#5) is shown with a dotted box delineating the la un- 
certainty in its photometric redshift. Radio-undetected SMGs are plotted at 
the most conservative low er limit on their redshift, given by the upper curve 
from lCarilli & Yunl feOOO) . Bottom: 450 fim/850 (im ratio as a redshift in- 
dicator. The curves are deriv ed from fits to the spectral e nergy distributions 
(SEDs) of z > 4 quasars JPriddev & McMahonll200ll solid) and SMGs 
iCoppin et al. 20071 dashed). 



Radio imaging at 1.4 GHz of the three fields, reaching rms lev- 
els of 17, 20 and 27/iJy beam -1 for SDSS J0756+4104, SDSS 
J1 030+0524 and SPS S J1044-0125, respectively, were published 
bv lPetric et all (T2003). As noted earlier, two sources in the SDSS 
J1030+0524 field have radio counterparts: identification of SMM 
J103026.29+052524.9 with a radio emitter implies z ~ 1, and 
SMM J103025.41+052431.7 is coincident with the core of a FRlI 
radio galaxy. Of the two quasars with submm detections, SDSS 
J0756+4104 has a 1.4-GHz flux density of 65 ± 17 /Jy, implying 
a spectral index, q? 5 4 °gh z z = 0.93 ± 0.05; SDSS J1044-0125 has 
a radio upper limit, 3cr < 80/iJy, giving af$=$t > 0.85 (2a). 
Both are consistent with starbursts, given their respective redshifts 
(Fig. |2j. In Fig. |2] we also indicate the lo wer limits on a^ 5 °g ^ 
corresponding to SMGs not detected in the lPetric e t al. (2003) ra- 



dio maps: all are consistent with the sources lying at high redshift, 
the weakest lower limits implying z > 1.5 (assuming a generous 
scatter in Qi^ghz •) 

Fig. [2] also illustrates the 450 /im/850 /im flux density ra- 
tio as a redshift indicator assuming idealised, isothermal, grey- 
body dust spectra: one SEP determined from z > 4 quasars by 
IPriddev & McMahonl ( 1200 ll) with T d = 42 K, /3 = 1.9; another 
SEP determined from blank-field SMGs in the SCUBA HAlf Pe- 
ree Extragalactic S urvey (SHAPES) with T d = 35 K, j3 = 1.5 
Coppin et al. 2007). A low 450 /im/850 /im flux density ratio im- 
plies a high redshift since the observed 450-/im emission must have 
originated at the peak, or possibly on the Wien side, of the SEP. 
Thus it is plausible that any 850-/im sources undetected in deep 
450-/im data lie at high redshift, the most conservative constraints 
giving z > 2. Conversely, a 450-/im emitter with no corresponding 
850-/im detection can be considered to have a weak upper limit on 
its redshift. 

Ideally, one would search for spectral features to establish 
the redshifts of the companion sources unambiguously. The rest- 
frame mid-IR, FIR and submm wavebands offer ionic/atomic fine- 
structure lines as well as the rotational transitions of molecular 
CO and HCN. Unfortunately, such observations are beyond the 
reach of current ground- and space-based instrumentation, though 
this will change with the advent of the Atacama Large Millime- 
tre Array (ALMA) and future actively cooled space missions such 
as the Spac e Infrared Telesc ope for Cosmology and Astrophysics 
(SPICA - iNakagawd l2004h or the FIR Interferometer (FIRI - 
Hel mich & Ivisonll2007t) .~ 



3.2 Number counts 

Earlier submm/mm imaging surveys have shown evidence for over- 
densities of SMGs in the vicinities of AGN at high redshift: in 

I 1 1 1 

fields around radio galaxies at z ~ 4 (e .g. Ivison et al. 2000a; 
IStevens et al.ll2003l : |Pe Breuck et ai]|2004|; iGreve et al. | l2007l) and 
in the field of an absorbed QSO at z ~ 2 dStevens et al.ll2004l) . Are 
our z > 5 QSO fields consistent with these findings? 

Fig. [H compares the cumulative number counts observed 
in o ur three z > 5 quasar fields ("z5Q") with (a) counts 
from lStevens et alj|2003l ("HzRG"); (b) bl ank-field sub mm survey 
counts (a fit to the SHAPES counts by ICoppin et alj|2006t): (c) 
counts from submm im aging surveys of clusters at z ~ 1 dBestl 
120021 ; IWebb et alj|2005l) . The raw number counts have been cor- 
rected for effective survey area as a function of limiting flux den- 
sity. However, we have not corrected for incompleteness or for flux 
boosting. At the faintest levels, or in noisy regions of the map where 
sources lie near the flux limit, these effects are significant, but they 
are difficult to correct in maps of this size and involve assump- 
tions about the number count distribution, e.g. using the blank-field 
counts with which we are attempting to compare. 

Nevertheless, assuming Poisson errors there is a significant ex- 
cess over the blank-field counts. The seven radio galaxy fields and 
the three z5Q fields each contain a total of 12 850-/im companion 
sources. The z5Q counts (and those of the z ~ 1 clusters) are con- 
sistent with the radio galaxy field counts (with the caveat that the 
z5Q sample spans a more limited range in flux density, due to shal- 
lower limit and smaller total area); both contain more SMGs than 
blank fields by a factor 4-5 across a range of 850-/im flux densities, 
including the ^6-mJy regime where the effects of flux boosting are 
minimal and overall sample reliability is excellent. 

In a blank-field survey of equivalent area and depth to all three 
of our fields combined, one would expect (based upon model fits to 
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Figure 3. Cumulative number counts of SMGs, in fields centred on (i) lumi- 
nous high-redshift AGN (filled symbols): z > 5 optically selected quasars 
(stars; this work); z ~ 4 radio galaxies (squares; Stevens et al. 2003); and 
(ii) z ~ 1 clusters (unfilled symbols - jBestil2002|;IWebb et alj|2005l) . The 
line shows a functional fit to the blank-field num ber counts of SMGs in 
the SHADES submm survey iCoppin et alj|2003) . Points have been cor- 
rected for effective area as a function of survey depth, but not for flux boost- 
ing, confusion or incompleteness. Even so, the excess over the blank-field 
counts is striking. 



the SHADES counts from ICoppinetal.ll2006h to detect ~3, ~2 
and ~1 sources brighter than Sgsoj™ = 4.1, 5.4 and 7.3 mJy, re- 
spectively. It is thus likely that some of the 12 z5Q companion 
sources are foreground contaminants, as we have seen already in 
£13.11 Others, as we have also seen, are consistent with SMGs at 
high redshift, based on their radio/submm and 450 /im/850 fim flux 
density ratios. We conclude that the maps probably contain SMGs 
genuinely associated with the dark matter haloes inhabited by our 
target quasars, leading to the statistical overdensities observed. 

As well as the potential for foreground objects to contaminate 
the sample directly via their submm emission, we should consider 
whether the observed exc ess of SMGs may be d ue to their gravita- 
tional lensing effects (e.g. lChapman et alj2002h. Depend ing on the 
slope of the luminosity function. IWvithe & Loebl d2002h estimate 
that ~ — ! of z ~ 6 quasars could be lensed. Targeting two 
submm-bright quasar hosts may have introduced an additional bias 
toward lensed fields, while the possibility that the submm-faint host 
is lensed is enhanced by the presence of a foreground radio galaxy 
(and possibly a cluster associated with that radio galaxy). The de- 
gree to which number counts are boosted by lensing depends on a 
competition between (i) the abundance of faint sources available to 
be boosted by the gravitational magnification and (ii) the stretch- 
ing of area in the source plane. Since the submm source counts are 
steep, (i) is a strong effect. Deep optical images of z > 5 quasars, 
including SDSS J1030+0524, ha ve failed to reve al any morpholog- 
ical signatures of strong lensing dFan et al.f2 003). Neither are line- 
of-sight galaxies, that could provide magnifications greater than 
~ 1.1, in evidence dWillottet al]2005l) . For now, therefore, we con- 
sider the observed overdensity to be due to high-redshift SMGs in 
the vicinity of the signpost quasars, rather than lensing or unusual 
foreground activity. 



3.3 Star-formation rates and dust masses at z > 5 

We can estimate SFRs for SMGs from their 850-^tm flux den- 
sities, S'ssoMm, assuming that they lie at the redshifts of the 



quasars. The range of FIR luminosity of a source at the aver- 
age redshift of the three quasars (z = 5.7) is Lfir = 0.7^1.4 
/mJy) x 1O 12 L , assuming thermal SEDs with Td = 
35-42 K and [3 = 1.5-1.9. If Lfir is powered by reprocessed stel- 
lar light, this translates into an instantaneous formation rate of mas- 
sive stars of * x (Lfir/1O 1o L0)Mq yr _1 . Here, ^ depends upon 
factors such as the stellar mass function and the efficiency with 
which starlight is reprocessed by dust; we assume <& ~ 1 giving 
a SFR range of 70— >140 x (5 , 85o f im/mJy) M Q yr _1 . Summing up 
the contributions of the companion SMGs (but not of the quasars, 
to avoid AGN contamination, even though it is likely that a sub- 
stantial fraction of their submm luminosity is due to star forma- 
tion) gives a total SFR per field, averaged over the three maps, of 
~ /z>sx 2000^4000 M yr -1 , where f z>5 represents the frac- 
tion of sources at the quasar redshifts. This star formation is appar- 
ently taking place in a region of volume 1 Mpc 3 . 

Adopting a d ust opacity Ki25 Mm = 30cm 2 g _1 (see 
IPriddev et al.l2003bl for an explanation), the dust mass correspond- 
ing to 5850/jm is M d = 0.5-^1.0 x (SssoMm/mJy) x 10 s Mq, giv- 
ing a total dust mass per field of w f z> $ x 1.6 — >2.8 x 10 9 M Q . In 
this case it is valid to include the quasar contribution, regardless of 
whether their submm luminosity is powered by starburst or AGN. 
The result is ~ 1.9 — >3.3 x 10 9 M ). If the dust-to-gas ratio in these 
objects is similar to that of the Milky Way, their total gas mass is 
~ (SgsoMm/mJy) x 1O 1O M . 

Finally, we consider the possibility that some of the submm 
companions are powered by buried AGN rather than obscured star 
formation. This is not implausible given evolutionary scenarios 
for SMGs involving co-evolution between black holes and mas- 
sive spheroids. Assuming, in the most extreme case, that Lfir 
represents the bolometric luminosity of the AGN, and assuming 
Eddington-limited accretion, the brightest SMGs (~10mJy) would 
have black holes of mass « 4 x 10 8 Mq . Their hard X-ray fluxes 
(2-10keV) would beRi4x 10~ 15 erg cm -2 s" 1 . For the general 
SMG population, or at least the radio-id entified subset, deep X-ray 
observations (e.g. Alexander et al. 2005) show that their AGN are 
weak compared with their bolometric luminosity. We have no rea- 
son to believe the SMGs in our high-redshift quasar fields should 
be any different; indeed, their S MBHs are likely somewhat less de- 
veloped than those explored bv lAlexander et all at z ~ 2.2. 



3.4 Resolved emission and radial profiles 

To measure source sizes requires knowledge of the PSF, determined 
to be 13.4 arcsec FWHM and near-circular by fitting a 2-dimensional 
(2-D) Gaussian to a beam map of 3C 345. We have used these 
parameters to determine the size of the most significant (> 6a) 
sources, again using 2-D Gaussian fits within the AIPS software en- 
vironment. 

SDSS J0756+4104 appears to be resolved along an axis with 
PA 69°. Along the orthogonal axis, the quasar's submm emission is 
point-like. Deconvolving the beam from the best-bet 2-D Gaussian 
fit suggests an emission region of size (16.0 ± 1.5) arcsec x (0.0 ± 
1.5) arcsec, though the apparent morphology could be mimicked by 
two or more well-separated, compact sources. 

At first sight SDSS 11044—0125 also appears resolved, al- 
beit with less certainty than for the ~14<r detection of SDSS 
J0756+4104. However, the best-fit 2-D Gaussian has dimensions 
similar to those of the beam: (14.4 ± 2.5) arcsec x (12.5 ± 
2.5) arcsec at PA 125°, so we conclude that there is no compelling 
evidence that the emission from SDSS J1044— 0125 has been re- 
solved. 



Submm environments of z > 5 quasars 7 



The angular scale is ~6 kpc arcsec at z ~ 5 — 6, so the phys- 
ical scale of the emission from SDSS J0756+4104 corresponds to 
~ 100 kpc on the plane of the sky. This suggests, perhaps, that we 
are witnessing a colossal merger. Lensing might provide an alterna- 
tive explanation, though there is scant evidence for this from optical 
images. 

Evidence of source structure deviating from that of a point 
source has been noted in previous submm observations of high - 
redshift galaxies (e.g. Ilvison et"al] [2000a; IStevens et af] 120031) . 
The mm continuum and CO line emission from the z = 4.7 
qu asar, BR 1202—0725 , was clearly separated into two components 
by lOmont etall dl996h . albeit with a smaller angular separation 
(4 arcsec) than single-dish submm imaging is capable of resolv- 
ing. Several high-reds hift radio galaxies have a l so been resolved 
in CO line emission ( Papadopoulos et al. 2000; De Breuck et al.l 
l2005h on arcsec scales using the IRAM Plateau de Bure interfer- 
ometer (PdBI). 

3.5 Comparison with model predictions 

The linear Press-Schechter approximation can be used to demon- 
strat e analytically the c oncept of bias in a hierarchical cosmology 
(e.g. lMo& White! l2002h . However, numerical simulations of the 
evolution of CDM, which can track the collapse of fluctuations on 
large scales into the non-linear regime, reveal a more complex pic- 
ture, the dark matter exhibiting a rich spatial structure (filaments, 

clusters: the so-called "cosmic web"). 

In the "Millenium Simulation", Sprin gel et alj d2005h specifi- 
cally identify z — 6 quasars, selecting them as the objects with the 
most massive dark matter haloes and/or the largest stellar mass. 
These objects have halo masses ~ 1O 12 ' 5 M , SFRs of several 
100 Mq yr _1 and evolve into the central dominant galaxies of rich 
clusters at z = 0. At z = 6 they are surrounded by numerous star- 
fo rming galaxies and l ie on prominent dark matter filaments (Fig. 3 
of lSpringel et alj2005h . The knots strung out on these filaments are 
reminiscent of the submm companions to z > 5 quasars - espe- 
cially, for example, the neighbouring pair in the SDSS J0756+4104 
image (Fig. [TJ, and the possible binarity of the quasar host itself. 
However, such claims are premature: the angular scale over which 
one would confidently expect to see large-scale structure is some- 
what greater than that enclosed within the field-of-view of SCUBA. 
The present images correspond only to the innermost ~ 1 Mpc of 
such structures. 

A tentative evolutionary scheme in whic h to interpret our ob- 
servations is the "anti-hierarchical" model of Granat o et al. 
2006). This scheme not only incorporates the effects of feedback 
from the AGN on star formation, but explicitly takes account of 
dust. A counter-intuitive feature of this model is that the evolution 
progresses more rapidly for the most massive objects. SMGs are en- 
visaged as massive spheroids undergoing a major episode of dust- 
enshrouded star formation, containing small (but gradually accret- 
ing) black holes in their cores. Eventually the black hole becomes 
sufficiently massive to power a quasar, thereby terminating star for- 
mation via jets in radio-loud systems or via accretion-driven winds 
in radio-quiet objects. The system then evolves passively as a mas- 
sive elliptical galaxy to the present day. Submm-luminous high- 
redshift quasars presumably corres pond to the late stages of this 
transi tion between SMG and QSO jPage et al. 2004; Stevens et al.l 
l2005h . If the QSO represents the most massive collapsed object in 
its field (by design we have selected extremely luminous objects), 
it will have a head start over the SMGs in its field. As suggested 
earlier, it may not be surprising that SDSS J1030+0524 is "submm 



dead", despite residing in the most active field: it may have evolved 
more rapidly than its companions and expended its fuel. 



4 SUMMARY 

Submm images at 450 and 850 fim of the fields of optically lu- 
minous z > 5 quasars reveal an excess of submm emitters rel- 
ative to expected counts of SMGs in blank fields. Although flux 
ratios show that some of the companion SMGs are undoubtedly 
foreground objects, the submm counts suggest that many of these 
companions probably lie in the dark matter haloes inhabited by our 
target quasars 

Two of the quasars are luminous submm emitters in their own 
right, suggesting that their hosts are intensely star-forming galax- 
ies. Furthermore, one shows tentative signs of extended submm 
emission, an indication either that the host galaxy is resolved (on 
a ~100-kpc scale) or, more likely, is caught during a merger. This 
finding reinforces, and extends to higher redshift, the conclusions 
of previous submm/mm imaging surveys of high-redshift AGN. 

The observed overdensities of luminous, star-forming galax- 
ies are consistent with the idea that luminous AGN reside in the 
most massive dark matter haloes at any epoch and, as such, pinpoint 
highly biased regions in the high-redshift Universe which eventu- 
ally merge to become rich clusters in the present day. Any more 
detailed interpretation is hindered by the limitations of present in- 
strumentation but forthcoming FIR and submm facilities will en- 
able significant progress to be made. 

• At current sensitivity levels, set by confusion at ~2mJy at 
850 /im, we are sensitive only to the very brightest starbursts - sev- 
eral xIOOMq yr _1 of star formation, not at all typical of known 
z > 5 field galaxies. However in the future, facilities such as 
ALMA will provide a powerful means of locating more typical 
star-forming galaxies at high redshift. It is important to charac- 
terise the feasibility of such a project by learning what we can about 
the role of dust in the formation of the earliest galaxies. The pres- 
ence of dust at these redshifts may have significant implications for 
reionization, as well as for Lyman-a line searches for high-redshift 
galaxies - the escape into the intergalactic medium of UV photons 
being hindered and thermally degraded by dust grains. 

• The large beam size of the JCMT at 850/iin precludes scrutiny 
of the sources on fine spatial scales, yet provides a tantalising 
glimpse of partially resolved or interacting galaxies. The enhanced 
resolution of submm interferometers such as IRAM PdBI and, ulti- 
mately, ALMA, will allow for more intricate study of the morphol- 
ogy and dynamics of high-redshift starbursts. 

• The high are al mapping speed offered by SCUBA-2 
jHolland et al]|2006l) opens up the possibility that we will be able 
to trace the distribution of star-forming galaxies on wide scales sur- 
rounding luminous, high-redshift AGN, obtaining surer statistics, 
perhaps even tracing the dark matter filaments and thus directly 
testing some of the predictions of numerical simulations of CDM 
models. 

• Wide-field spectral imaging at mid-IR-mm wavelengths of- 
fers the best prospects with which to determine the redshifts of 
SMGs in the quasar fields unambiguously. Whilst such observa- 
tions are beyond the capabilities of current facilities, the sensitivi- 
ties of ALMA and future FIR space missions will be sufficient to 
pinpoint dusty, star-forming galaxies out to z ~ 6. 

The prospects implied by these discoveries for the next gener- 
ation of submm instruments are therefore extremely promising. 
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